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Breaking of the number-of-constituent-quark scaling for identified-particle elliptic flow 
as a signal of phase change in low-energy data taken at the BNL Relativistic Heavy 

Ion Colhder (RHIC) 
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We argue that measurements of identified-particle elliptic flow in a wide energy range could 
shed light on the possible phase change in high-energy heavy ion collisions at the BNL Relativistic 
Heavy Ion Collider (RIfIC). When the hadronization process is dominated by quark coalescence, 
the number-of-constituent-quark (NCQ) scaling for the identified-particle elliptic flow can serve as a 
probe for studying the strong interacting partonic matter. In the upcoming RHIC low-energy runs, 
the NCQ scaling behavior may be broken because of the change of the effective degrees of freedom 
of the hot dense matter, which corresponds to the transition from the dominant partonic phase to 
the dominant hadronic phase. A multiphase transport model is used to present the dependence of 
NCQ scaling behavior on the different hadronization mechanisms. 
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The main purpose of the heavy ion program at Rel- 
ativistic Heavy Ion Collider (RHIC) in Brookhavcn Na- 
tional Laboratory (BNL) is to study the properties of the 
high density matter produced in the heavy ion collisions, 
in particularly, whether it undergoes a transition from 
hadronic phase to quark gluon plasma (QGP) phase. Re- 
cently, many experimental results demonstrated that a 
hot and dense partonic matter has been indeed formed 
in ultra-relativistic energy heavy ion collisions at ^/sJjn 
= 200 GeV [H . In order to study this new form of mat- 
ter, probes are essential to gain information from the 
earliest stage of the collisions. Measurements of the col- 
lective motion, especially the elliptic flow V2 of identified 
particles produced in heavy ion collisions have long been 
suggested as a valuable tool to study the nature of the 
constituents and the equation of state of the system in 
the early stage of the reaction Many significant re- 
sults in low transverse momentum range for elliptic fiow 
V2 of final state particles have been obtained from RHIC 
experiments d, H, [E H 13 which consist well with the 
predictions from ideal hydrodynamics Q. More inter- 
estingly, a number of constituent quark (NCQ) scaling 
behavior has been observed in the intermediate pr range 
(1.5 GeV/c < pt < 5 GeV/c), which can be well repro- 
duced by parton coalescence and recombination model 
calculations [l^, [111, [l^- Such scaling indicates that 
the collective elliptic flow has been developed during the 
partonic stage and the effective constituent quark degree 
of freedom plays an important role in the hadronization 
process. 

Simultaneously, some experimental and theoretical de- 
velopments have suggested that important discoveries 
are possible at lower collision energies where physics 
base may be completely different from the pictures at 
higher energies. Specifically, calculations from lattice 
QCD [13, 14] predict a transition or fast cross-over be- 
tween the QGP state and the hadronic matter at Tc ~ 
150'--^180 MeV with vanishing baryon density. In ad- 



dition, several observables in central Pb-|-Pb collisions 
from SPS experiments show qualitative changes in the 
energy dependence [lH|, e.g., the full phase space ratios 
{K+)/{tt+) and ((A) -I- {K + 1{))/{tt) showed a turnover 
and a decrease around y^sjvw ~ 7.6 GeV after a steep 
increase at lower energies. 

To gain the further information about the existence 
of the possible phase transition between hadronic matter 
and partonic matter, the future RHIC physics program 
includes the Au-|-Au energy scan extending to low colli- 
sion energies. Based on this, we shall focus our studies 
on the NCQ-scaling of the identified particles elliptic fiow 
in a wide energy range. The turning on/off of the NCQ- 
scaling behavior may indicate the onset of deconfinement: 
the scaling will be retained in the partonic phase at high 
energy, while it may be broken in lower energy when the 
system is dominated by hadronic interaction. 

In the non-central collisions, the initial asymmetries in 
the geometry of the system can lead to the anisotropics of 
the particle momentum distributions. Since the spatial 
asymmetries decrease rapidly with time, anisotropic fiow 
can develop only in the first few fm/c. In that way, the 
properties of the hot dense matter formed during the 
initial stage of heavy ion collisions can be learned by 
measuring the anisotropic fiow. The anisotropic flow is 
defined as the nth Fourier coefficient w„ of the particle 
distributions in emission azimuthal angle with respect to 
the reaction plane [l^ , which can be written as 



— — cx 1 + 2 > Vn cos(nA0j, 
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where A0 denotes the angle between the transverse mo- 
mentum of the particle and the reaction plane. The 
second Fourier coefficient V2 represents the elliptic fiow 
which characterizes the eccentricity of the particle distri- 
butions in momentum space. At a given rapidity window 
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FIG. 1: (Color online) Experimental results of elliptic flow as a 
function of transverse momentum for tt, Kg, p, A, </>, S, and O with 
\ri\ <1.0 in 0-80% Au+Au collisions at ^/sJTn = 200 GeV from 
RHIC. The results are taken from STAR Experiments HUHIll. 



FIG. 2: (Color online) Number of constituent quark (riq) scaled 
V2/nq versus scaled (my — m)/nq. Dashed blue line is the result 
of fit with function (3). 



the second coefficient is 

i;2 = (cos(2A0)> = (4^). (2) 

The measured elliptic flow V2 as a function of transverse 
momentum px from the minimum bias Au+Au collisions 
at y/sNN = 200 GeV for tt, Kg, p, A, S, cj), Q are pre- 
sented in FIG.[TJ The results show that elliptic flow V2 of 
all the particles increase with the transverse momentum 
while the elliptic flow V2 becomes saturated in the in- 
termediate transverse momentum region. Moreover, the 
baryons saturate at px > 5 GeV/c with V2 ~ 0.2, while 
the mesons saturation starts earlier at lower values of 
V2- In addition, from FIG. [T]we can find that although 
the multi-strange baryons S and tend to undergo less 
re-scatterings in the hadronic stage, their V2 values are 
as high as others hadrons at given px range, this means 
that the collectivity should be developed at the partonic 
stage. 

More intriguing phenomenon is that after scaling both 
values of V2 and px by the number of the constituent 
quarks of the corresponding hadron, all particles can fall 
onto one single curve. However, the pion elliptic flow 
data are somewhat not following the scaling, which may 
be caused by the large contribution of pion yield from res- 
onance decays [13) Ell- FIG. [5] shows elliptic flow V2 as 
a function of transverse kinetic energy {mx — m) for the 
identified particles {nix — \J + Px is the transverse 
mass while m is the rest mass of the particle) , where V2 
and (jnx — m) have been scaled by the number of con- 
stituent quarks {uq). The measured elliptic flow values 
can be fitted with the equation 18*] given as 

^■"^ = l + exp(-(x-fo)/c) " 
where a, 6, c and d are the parameters fixed from the 



fit. This universal curve represents the momentum space 
anisotropy of constituent quarks prior to hadron for- 
mation. However, this simple scaling neglects possible 
higher harmonics and possible differences between light 
and heavier quark flow. 

Some coalescence or recombination models d, [13, [ii[, 
can successfully describe the hadron production in 
the intermediate px region, where the NCQ-scaling be- 
havior has been observed. According to these models, 
the hadronization is dominated by coalescence of quarks, 
and the essential degrees of freedom seem to be effective 
constituent quarks that have developed a collective ellip- 
tic flow during the partonic evolution. There are some 
common features in the intermediate px region below 5 
GeV/c [l^: (i) The production probability for a baryon 
or meson is proportional to the product of local parton 
densities for the constituent quarks; (ii) Baryons with 
transverse momentum px are mainly formed from quarks 
with transverse momentum ~ Pt/-^, whereas mesons at 
Px are mainly produced from partons with transverse 
momentum Pt/'^- Then the meson elliptic flow V2,m 
and baryon elliptic flow V2,b can be given by those of 
partons via (20,] 

Px 

V2,MiPT) ~ 2v2,q{ — ), (4) 
V2,b{Pt) ~ 3w2,g(^)- (5) 

Consequently, the collectivity of the constituent quarks 
become the collectivity of the hadrons via quark coales- 
cence during the hadronization. Therefore, by scaling the 
observed V2 signal and the transverse momentum with 
the number of constituent quarks Uq, one obtains the 
underlying quark flow. 

Many experimental results [l5| have shown anomalous 
dependence on the collision energy in the SPS energy 
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FIG. 3: (Color online) (Left panel) Number of constituent quark (riq) scaled V2/nq versus scaled (m-r — m)/nq at Au+Au y'sjvjv= 12.3 
GeV from the AMPT model with string melting. Dashed blue line is the result of fit with the function (3). (Right panel) Number of 
constituent quark (riq) scaled V2/nq versus scaled (mj' — m)/nq at Au+Au -y/sjvjv= 12-3 GeV from the AMPT model without string 
melting. 



range, especially the sharp changes around ^/snn ~ 7.6 
GeV. All those phenomena suggest that a common un- 
derlying physics process is responsible for these changes. 
Ref. [2l|, [24I argue that this anomaly is probably caused 
by the modification of the equation of state in the transi- 
tion region between confined and deconfined matter. The 
low energy runs at RHIC offer us opportunities to study 
the possible phase transition, thus the NCQ-scaling be- 
havior for the identified particles may be considered as 
a new signal of the onset of deconfinement located in 
the low energy domain. We suppose that at lower col- 
lision energies, there would be no phase transition from 
hadronic matter to partonic matter. In this case, the 
hadrons are not produced from the coalescence of decon- 
fined partons, and the collective flow does not stem from 
partonic stage, therefore the NCQ-scaling behavior for 
identified particles will be broken. 

In the following, we will use a multi-phase trans- 
port (AMPT) model ^ to investigate the NCQ- 
scaling behavior of the identified particles under dif- 
ferent hadronization mechanisms in Au-|-Au collisions. 
The AMPT model has an extensively agreement with 
the RHIC data and AGS data [H, [3, l25, ^fi, ^ 113. 
For example, by using parton-scattering cross sections 
of 6-10 mb, it was able to reproduce both the centrality 
and transverse momentum (below 2 GeV/c) dependence 
of the elliptic fiow and pion interferometry measured in 
Au-|-Au collisions at ^/snn = 130 GeV. It can explain the 
measured px dependence of both V2 and V4 of midrapidity 
charged hadrons and V2 for meson as well as di-hadron 
correlation in the same collisions at y^s^N = 200 GeV. It 



was also demonstrated that there exists the NCQ-scaling 
of V2 not only for tt, k and p but also for multi-strange 
hyperons in the AMPT model with the string melting 
scenario [H, . The AMPT model is a hybrid model 
which contains many processes of Monte Carlo simula- 
tion. There are four main components in the model: 
initial conditions, partonic interactions, conversion from 
partonic matter into hadronic matter, and hadronic in- 
teractions in collision evolution. First the hybrid model 
uses the minijet partons from the hard processes and the 
strings from the soft processes in the HIJING model as 
initial phase j29| . then the dynamical evolution of par- 
tons are modeled by the ZPC [30[ parton cascade model, 
which calculates two-body parton scatterings using cross 
sections from pQCD with screening masses. The tran- 
sition from the partonic matter to the hadronic matter 
is based on the Lund string fragmentation model [3]| . 
The final-state hadronic scatterings are modeled by the 
ART model 32]. In the default AMPT model, minijets 
coexist with the remaining part of their parent nucleons, 
and together they form new excited strings, then the re- 
sulting strings fragment into hadrons according to the 
Lund string fragmentation. In the AMPT model with 
the string melting scenario, these strings are converted 
to soft partons, and their hadronization is based on a 
naive quark coalescence model. 

We calculated the elliptic flow V2 with different ver- 
sions of AMPT model in Au-|-Au collisions at ^/snn = 
12.3 GeV, which assumed to be running at RHIC in the 
future. Based on the momentum of identified particles, 
we extract the elliptic flow V2 with formula (2). Accord- 
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ing to the discussion above, when there is a phase tran- 
sition from hadron matter to quark gluon plasma, the 
essential degrees of freedom at the hadronization seem to 
be effective constituent quarks, then we use the AMPT 
model with string melting to calculate the elliptic flow. 
Left panel of FIG. [3] shows the elliptic flow V2 for iden- 
tified particles calculated with string melting scenario, 
it seems that the elliptic flow V2 of the identifled parti- 
cles have a nice NCQ-scaling behavior due to the coales- 
cence of partons in the hadronization. On the other hand, 
if the degree of freedom of the system is dominated by 
the hadronic interaction, we will use the default AMPT 
model to calculate the elliptic flow, right panel of FIG. [3] 
shows the elliptic flow V2 of identifled particles without 
string melting, the NCQ-scaling of elliptic flow V2 is bro- 
ken since hadrons are directly from the fragmentation of 
excited strings. In comparison with the string melting 
case, the V2 values of identifled particles show significant 
decreasing in the case without the string melting, which 
refiects that partonic stage makes an important contri- 
bution on the development of sizable elliptic fiow in the 
early stage. 



In summary, we proposed to use the NCQ-scaling for 
identified particles elliptic fiow as a unique observable to 
probe the effective degree of freedom of the system cre- 
ated in heavy ion collisions at RHIC. We argued that the 
presence or absence of the NCQ-scaling behavior might 
indicate the turn on or turn off of the partonic degree 
of freedom of the system. The AMPT simulation has 
further confirmed our argument: the NCQ-scaling is ob- 
served in the string melting scenario where full partonic 
evolution has been included while it is broken in the de- 
fault version with dominant hadronic interaction only. 
The upcoming RHIC low energy measurement will help 
to disentangle the diff'erent physics scenario as discussed. 
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